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Introduction 


As a result of the ongoing exploration of Mars and the several unmanned and possibly 
manned missions planned for the near future, increased attention has been given to the use of the 
natural resources of the planet for rocket propellant production and energy generation. Since the 
atmosphere of Mars consists of approximately 95% carbon dioxide (CO 2 ), this gas is the resource 
of choice to be employed for these purposes. Since many metals burn vigorously with CO 2 , these 
may be used as an energy source or as propellants for a research vehicle on the surface of Mars. 
Shafirovich and Goldshleger 1 conducted experiments with spherical particles up to 2.5 mm in 
diameter and found that the burning process was controlled by diffusion and that the particles 
exhibited pulsating combustion due to superheating of the Mg vapor trapped inside a protective 
oxide shell. They also proposed a reaction mechanism based on the gas-phase reaction. Mg + 
CO 2 —> MgO + CO and the heterogeneous reaction Mg + CO — » MgO + C occurring on the 
sample surface. 

In all the above studies with large Mg particles, the burning process is invariably 
influenced by strong convective currents that accelerate the combustion reaction and shorten the 
burning times. Although these currents are nearly absent in the burning of small particles, the high 
emissivity of the flames, rapid reaction, and small length scales make the gathering of any useful 
information on burning rates and flame structure very difficult. The goal of this investigation is to 
provide a detailed study of flame structure by taking advantage of large, free-floating spherical 
metal samples and their corresponding long burning times available in a weightless environment. 
The use of reduced gravity is essential to eliminate the intrusive buoyant flows that plague high- 
temperature metal reactions, to remove the destructive effect of gravity on the shape of molten 


metal samples, and to study the combustion behavior of metals in the presence of solid oxides 
undisturbed by natural convection. This work presents the most complete modeling of metal 
particle burning to date for Mg with CO2 and O2. 

Numerical Model 

In the case of metal combustion experiments conducted in microgravity, a spherical 
symmetry during burning and the absence of natural-convection buoyant plumes are nearly 
achieved. The present experimental effort is complemented by the development of a numerical 
simulator, DROPLET, to model the spherically symmetric, quasi-steady burning of a metal fuel 
droplet with full gas-phase chemical kinetics and transport properties calculations. Due to the 
lack of kinetic and thermophysical data for the heterogeneous reactions of Mg and the absence of 
a complete condensation mechanism of MgO, the present model does not take into account 
surface reactions, radiation heat transfer, or complete condensation processes. Only vaporization 
of Mg from the sample surface is considered (with no heterogeneous surface reactions) and MgO 
condensation is treated with an Arrhenius-type rate expression for the molar production rate of 
the liquid, MgO(l) and solid oxide, MgO(s). The condensed products, MgO(l) and MgO(s), are 
treated like gas-phase species with a low diffusivity. 

There are several numerical models already developed for the combustion of aluminum 
(Al) particles 2 " 5 that take into account some of the above processes not included in the present 
simulation. Unfortunately, these models can not be used in the case of Mg combustion since 
knowledge of the kinetic and thermophysical properties of Mg and its oxide is not as extensive as 
for Al. Thus, the present model was developed to help in the interpretation of the experimental 


results and to serve as a preliminary effort for further development of a more complete and much 
needed model of Mg combustion. 

Before attempting a Mg-C0 2 flame simulation, preliminary simulations of Mg-0 2 and Mg- 
air flames are performed since the reaction mechanisms for these systems are better understood. 
Equilibrium calculations by the NASA-Lewis chemical equilibrium code 6 are performed for all the 
above metal-oxidizer combinations to identify the most important species and to obtain the 
equilibrium temperature. 

Mathematical Formulation 

The one-dimensional conservation equations describing the quasi-steady burning of a spherical 
droplet in a quiescent gaseous atmosphere are: 
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The independent variable in these equations is the radial coordinate, r. The mass density, 
p, is calculated from the ideal-gas equation. The model assumes a uniform droplet temperature 
and the pressure is considered constant throughout the domain. 

The multicomponent transport formulations are used for calculation of all transport 
properties. The species diffusion velocities are then expressed as: 
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Boundary conditions 


In order to solve the governing equations described above the following boundary 
conditions are imposed: 
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Equation (5) is the species mass balance for the species at the surface implying that the 
chemical species flux at the surface is the same as the amount produced or destroyed by the 
surface chemical reaction. Equation (6) is the surface energy balance and Eq. (7) is the Clausius- 
Clapeyron liquid-vapor equilibrium relationship applied to the fuel. 

In solving the above equations it is assumed that there is no reaction at the surface, so the 
right-hand side is zero in Eq. (5). Equation (1) can be simplified as: 


r 2 pu r = M f 


( 8 ) 


where the mass flux, M h is a constant. Thus, M f is conserved throughout the domain. Use of this 
conserved scalar simplifies other governing equations. Hence, there is no need to solve for the 
velocity field. 

Solution Method 

Finite-difference approximations for non-uniform mesh spacing are used to solve the 
resulting system of partial differential equations by forming corresponding algebraic equations. 
Diffusive terms use central differences with truncation error that is second order in mesh spacing. 
For better convergence, convective terms use upwind differencing with truncation error that is 
first order in the mesh spacing. Use of Eq. (8) suggests that the resulting matrix structure is non- 
symmetric. Since TWOPNT 7 , the boundary value problem solver used by this simulator, needs a 
block tridiagonal matrix structure, this equation is copied at each grid point. Also, no boundary 
condition is required for the mass flux far from the particle. The basic solution method is based 
on a damped modified Newton iteration. A time stepping algorithm is incorporated to help bring 
the successive iterations within the domain of convergence for Newton's method. 

DROPLET is the simulator written for this specific problem. It uses CHEMKIN 8 and the 
transport package for kinetics and transport calculations, respectively. The flowchart of the 
simulator is shown in Fig. 1. The input to CHEMKIN consists of the reaction mechanism and the 
thermodynamic data. Input to the transport package is the transport property data. The 
CHEMKIN interpreter and the transport fitting code write respective linking files that can be 
retrieved by the corresponding libraries, such as CKLIB and TRANLIB. It also uses the 
boundary value problem solver, TWOPNT, to solve the resulting system of algebraic equations. 
The input file to DROPLET takes input from the user that consists of the pressure, initial number 


of grid points, ambient (far field) conditions, and guessed values of surface and maximum 
temperatures. The output file written by DROPLET contains the initial guess and the final 
solution along with the intermediate solutions found on smaller grids. 

Numerical Model Predictions 

As explained in the Numerical Model section, numerical simulations are first conducted 
for the Mg-02 and Mg-air systems for which the reaction mechanisms are better understood. The 
reactions in the Mg-0 2 system are then used in the Mg-C0 2 simulation along with the reactions of 
Mg with C0 2 and CO, as well as all reactions of carbon-containing species with 0 2 and O. For all 
three metal-oxidizer systems mentioned above, equilibrium calculations are performed with the 
NASA-Lewis chemical equilibrium code 6 to identify the most important species and to obtain the 
equilibrium temperature. For each of the three systems, it is necessary to include all three phases 
of MgO (gas, liquid, and solid) to obtain a realistic equilibrium temperature as observed in 
experiments. Including the gas- and liquid-phase species, MgO and MgO(l), as the only oxide 
products results in unrealistically low temperatures and species concentrations due to the absence 
of the large heat release experienced from the complete gas-to-solid condensation process. 
Hence, both the gas-to-liquid and liquid-to-solid condensation mechanisms must be considered for 
an accurate description of the Mg combustion process. 

The Mg-0 2 flame structure 

In the case of Mg-0 2 , the equilibrium calculations are performed with 1 mole of gaseous 
Mg at 1366 K and 1 mole of 0 2 at 300 K as the initial conditions. The equilibrium temperature 


under adiabatic conditions is 3398 K and all three phases of MgO are present at this temperature. 
The numerical simulation of the Mg-02 flame structure uses the following reaction mechanism. 


Mg + 02 = MgO + O 

(Rl) 

Mg + 0 + M = MgO + M 

(R2) 

MgO = MgO(l) 

(R3) 

Mg 0(1) = MgO(s) 

(R4) 

20 + M = 0 2 + M 

(R5) 


Reactions (Rl) and (R2) and their corresponding rate constants are determined from 
previous studies 9 . The symbol M represents the presence of a third body. The molar production 
rate of MgO(l) in reaction (R3) is obtained from the Arrhenius- type rate expression 10 of the 
condensation rate for the formation of drops of the critical radius (condensation nuclei) arising in 
a unit volume during a unit time. Due to the lack of a similar rate expression for the licjuid-to- 
solid phase transition, a high value of the pre-exponential factor is used to simulate the 
condensation process of reaction (R4). This parameter is adjusted until a physically acceptable 
temperature profile is achieved. The condensed products are treated as gas-phase species with 
very low diffusivity. This is achieved by using a large value of the Lennard- Jones collision 
diameter in the transport property input. In the present simulations a Lennard-Jones collision 
diameter of 10 A is chosen because is of the order of magnitude of the critical radius for 
condensation (3 A) of MgO under the given conditions and small enough to avoid the treatment 
of a condensed phase. This value is used as a collision diameter for both MgO(l) and MgO(s). 
The rate constants for reactions (Rl) through (R5) are shown in Table 1. 

Figure 2 shows the flame structure of a 2-mm diameter Mg droplet burning in pure O 2 at 1 
atm pressure and 300K temperature. Magnesium vapor diffuses out from the particle surface 


towards the counterdiffusing 0 2 . The surface temperature is only a few degrees below the boiling 
point of Mg (1366K). The temperature profile rises with a very steep gradient near the surface to 
a short plateau region where a maximum temperature of 3220K is reached. After the complete 
consumption of the Mg vapor, the slope of the temperature profile starts decreasing at a faster 
rate. The predicted maximum temperature of 3220K is slightly lower than the maximum 
equilibrium temperature of 3398K, i. e. the vaporization-decomposition point of MgO under these 
conditions. This difference in maximum temperature values may be attributed to the lack of an 
accurate condensation model which may also be responsible for the absence of a large 
concentration of MgO(l) in the plateau region where the temperature exceeds the melting point of 
MgO (3105 K) and where the liquid oxide is the dominant species. The high pre-exponential 
factor used to simulate the liquid- to -so lid phase change is insensitive to temperature and thus 
unable to simulate the rise of MgO(l) in the temperature plateau region. Hence only a very small 
amount of MgO(l) is obtained in the present calculations. 

The gaseous oxide profile, MgO, shows a peak at the highest temperature and then rapidly 
decreases. As a result of the high flame temperature, a large amount of atomic oxygen, O, is also 
found in the plateau region of the temperature profile. The simulation also shows the diffusion of 
molecular oxygen, 0 2 , to the particle surface. Its significant concentration at the surface indicates 
that heterogeneous reactions may play an important role at the surface of the sample. 

As mentioned before, the liquid and solid oxides are treated as gas-phase species with low 
diffusivity by using a Lennard-Jones collision diameter of 10 A, a value of the same order of 
magnitude as the critical radius of condensation. To test the effectiveness of this approach, a 
collision diameter of 25 A was used under the same conditions as in the original simulation. 
Under this higher collision diameter, the burning rate decreases by a factor of two and hence the 


burning time increases by the same magnitude. Therefore, it appears that increasing the 
agglomeration of condensed products around the reaction zone results in increased hindrance to 
the transport of the gaseous oxidizer. However, since the estimate of burning times and hence the 
corresponding d 1 law are strongly dependent on the assumed value of the collision diameter, no 
direct comparison with experimental burning times can be done without the inclusion of complete 
condensation and agglomeration mechanisms. Nevertheless, as will be seen in the following 
sections, comparisons can be done between burning rates obtained by simulation for the Mg-0 2 , 
Mg-air, and Mg-C0 2 systems. 

The Mg- Air flame structure 

In the case of the Mg-air system, the equilibrium calculations are performed with 1 mole 
of gaseous Mg at 1366 K and 1 mole of air at 300 K as the initial conditions. The equilibrium 
temperature under adiabatic conditions is 3220 K with all three phases of MgO present at this 
temperature. For the numerical simulation of the Mg-air flame structure, nitrogen, N 2 , is added 
(79% by volume) to the oxidizer to simulate the burning of Mg in air. The GRI-Mech™ 2.11 
mechanism 11 is used to add the nitrogen chemistry to the Mg-0 2 mechanism. Other product 
species such as Mg 3 N 2 and MgN are neglected in the present analysis due to the lack of kinetics 
and transport property information. The condensation of MgO into MgO(l) and MgO(s) is 
treated in the same manner as in the Mg-0 2 case. The reaction mechanism and rate constants for 
the Mg-air system are shown in Table 2. 

Figure 3 shows the flame structure of a 2-mm diameter Mg droplet burning in air at 1 atm 
pressure and 300K temperature. As in the Mg-0 2 flame, the temperature profile rises sharply near 
the droplet surface. The maximum temperature in this case, however, is 3143 K, which is lower 


than in the Mg-0 2 system due to the presence of N 2 . This predicted maximum temperature of the 
Mg-air flame is also lower than the adiabatic flame temperature obtained in the equilibrium 
calculations (3220K) due to the lack of an accurate condensation model, which is also responsible 
for the absence of any significant concentration of MgO(l) in the temperature plateau region. The 
predicted burning rate is 1.5 times smaller than that obtained in the Mg-0 2 simulation. Also, an 
increase in the Lennard-Jones collision diameter to 25A decreases the burning rate by a factor of 
two, which is in excellent quantitative agreement with the trend obtained in the simulation of Mg 
burning in pure oxygen. 

The Mg-C0 2 flame structure 

For the Mg-C0 2 system, an equilibrium calculation performed with 1 mole of gaseous Mg 
at 1366 K and 1 mole of C0 2 at 1000 K gives an equilibrium temperature under adiabatic 
conditions of 3 174 K with all three phases of MgO present at this temperature. For the numerical 
simulation of the Mg-C0 2 flame structure, the reaction mechanism proposed by Shafirovich and 
Goldshleger 1 (described in the Introduction) is used for the reactions of Mg with C0 2 and CO: 

Mg + C0 2 = MgO + CO (R 11) 

Mg + CO = MgO + C (R12) 

The elementary steps and rate constants for reactions of carbon-containing species with 0 2 
and O are taken from the GRI-Mech™ 2.1 1 mechanism" and are shown in Table 3 along with all 
the Mg-0 2 and Mg-C0 2 reactions. Since the rate constants for reactions (R1 1) and (R12) are not 
well known, these constants are obtained for the present study from a sensitivity analysis and its 
comparison with experimental results. Species such as MgCOj, MgC, and Mg 2 C 3 , predicted by 


the equilibrium calculations, are neglected in the present analysis due to the lack of kinetics 
information. Condensation processes are treated in the same manner as in the Mg-O; and Mg-air 
cases. 

The flame structure of a 2-mm diameter Mg droplet burning in a pure CO 2 atmosphere at 
1 atm pressure is shown in Fig. 4. In this case the oxidizer ambient temperature is 1000 K since 
no convergence is obtained in the simulation for lower temperatures. This numerical result is in 
agreement with experimental studies 1 where the Mg sample is ignited by slowly heating the CO 2 
gas around it. In these tests, no burning is obtained for ambient temperatures below 985 K. 
Under the experimental conditions used in the present investigation, the burning of the Mg sample 
under CO 2 ambient temperatures as low as 300 K is possible because enough radiant energy is 
provided to the sample to induce its ignition and subsequent burning before the energy source is 
turned off. The maximum combustion temperature predicted by the simulation is 2645 K, which 
is lower than the equilibrium temperature of 3174 K for CO 2 at an ambient temperature of 1000 
K. The cause of this discrepancy in maximum temperature values follows the same reasoning as 
in the Mg-0 2 and Mg-air cases. 

The reaction zone observed in Fig. 4 is narrower than in the Mg-C >2 and Mg-air 
simulations, suggesting a slower molecular diffusion process in the case of C0 2 as the oxidizer. 
This result is in qualitative agreement with the experiments where longer burning times are 
observed in Mg-C0 2 flames. The predicted burning rate is slightly lower than the rate obtained in 


the Mg-air simulation. 


Conclusions 


This investigation studies the burning behavior and flame structure of Mg in a C0 2 
atmosphere to assess the feasibility of using metal-C0 2 reactions as an in situ resource utilization 
technology for rocket propulsion and energy generation on other planets. A one dimensional, 
spherically symmetric, quasi-steady model is used to simulate the burning behavior of Mg in 0 2 , 
air, and C0 2 . The flame structure predicted by the model shows a temperature profile with a 
short plateau region at the highest temperature of the flame. It is within this region that the 
vaporization-dissociation of MgO takes place acting as an enthalpy sink. Condensation of MgO 
to MgO(l) and MgO(s) follows thereafter as it diffuses away from the flame. The maximum 
temperature obtained from the model is lower than that predicted by the equilibrium calculations 
and the high concentration of MgO(l) expected in the temperature plateau region is absent from 
the simulation results. Both effects are due to the lack of accurate information on the rates of 
liquid-to-solid condensation of the metal oxide. 

In accordance with the experiments, the overall burning process appears to be diffusion 
controlled in all the cases. The burning rate and the maximum temperature in a C0 2 atmosphere 
are lower than in 0 2 or air. For Mg burning in C0 2 , an ambient temperature of 1000 K is 
necessary for convergence of the simulation, which agrees with experimental observations of 
samples ignited and burned by slow heating of the surrounding oxidizer. 

At the present time, this preliminary model does not give a complete picture of the metal- 
oxidizer reaction since not all of the relevant physical and chemical mechanisms are considered. 
Product condensation and agglomeration, radiation heat loss, and heterogeneous surface reactions 
are just some of the important mechanisms that should be included in the model to provide an 
accurate representation of the burning process. In addition, there is also a lack of information in 


the literature on reaction rates of the most basic elementary reactions. Nevertheless, although 

preliminary, the model provides a qualitative description of the flame structure that helps to 

visualize the complex phenomenon of metal combustion. 
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Figure 1 . Flowchart of the DROPLET numerical simulator. 
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Table 1. Reaction Mechanism for Mg-0 2 


Reaction 


A 

P 

e a 

Ref. 

Mg + 0 2 = MgO + 0 

(Rl) 

4.44E+12 

0.5 

30500.0 

9 

Mg + 0 + M = MgO + M 

(R2) 

1.90E+14 

0.5 

0.0 

9 

MgO = MgO(l) 

(R3) 

6.11E+11 

-2.0 

17268.6 

10 

MgO(l) = MgO(s) 

(R4) 

l.OOE+15 

0.0 

0.0 

— 

20 + M = 0 2 + M 

(R5) 

1.89E+13 

0.0 

-1788.0 

12 

The rate coefficients are in the form Kj ■ 

= A* T* 3 Exp(-E A /RT). Units are in moles, 

cubic centimeters, 

seconds, 

Kelvins, and calories/mole. 






Table 2. 

Reaction Mechanism for Mg-Air 



Reaction 


A 

P 

e a 

Ref. 

Mg + 0 2 = MgO + O 

(Rl) 

4.44E+12 

0.5 

30500.0 

9 

Mg + 0 + M = MgO + M 

(R2) 

1.90E+14 

0.5 

0.0 

9 

MgO = MgO(l) 

(R3) 

6.11E+11 

-2.0 

17268.6 

10 

MgO(l) = MgO(s) 

(R4) 

l.OOE+15 

0.0 

0.0 

— 

20 + M = 0 2 + M 

(R5) 

1.89E+13 

0.0 

-1788.0 

12 

N + NO = N 2 + 0 

(R6) 

3.50E+13 

0.0 

330.0 

11 

N + 0 2 = NO + 0 

(R7) 

2.65E+12 

0.0 

6400.0 

11 

n 2 o + 0 = n 2 + o 2 

(R8) 

1.40E+12 

0.0 

10810.0 

11 

N 2 0 + O = 2N0 

(R9) 

2.90E+13 

0.0 

23150.0 

11 

N0 2 + 0 = NO + 0 2 

(RIO) 

3.90E+12 

0.0 

-240.0 

11 

The rate coefficients are in the form Kj 

= A*T ** Exp(~E A /RT). Units are in moles, 

cubic centimeters, 

seconds, 


Kelvins, and calories/mole. 


Table 3. Reaction Mechanism for Mg-C02 


Reaction 


3 

P 

E a 

Ref. 

Mg + C0 2 = MgO + CO 

(Rll) 

2.00E+14 

0.5 

34847.0 

1 

Mg + CO = MgO + C 

(R12) 

2.00E+11 

0.5 

30000.0 

1 

Mg + 0 2 = MgO + O 

(Rl) 

4.44E+12 

0.0 

30500.0 

9 

Mg + 0 + M = MgO(s) + M 

(R2) 

1.90E+14 

0.0 

0.0 

9 

MgO = MgO(l) 

(R3) 

6.11E+11 

-2.0 

17268.6 

10 

MgO(l) = MgO(s) 

(R4) 

1.00E+15 

0.0 

0.0 

— 

20 + M = 0 2 + M 

(R5) 

1.89E+13 

0.0 

-1788.0 

12 

CO + 0 + M = C0 2 + M 

(R13) 

6.02E+14 

0.0 

3000.0 

11 

CO + 0 2 = C0 2 + O 

(R14) 

2.50E+12 

0.0 

47800.0 

11 

C + 0 2 = CO + O 

(R15) 

5.80E+13 

0.0 

576.0 

11 


The rate coefficients are in the form K } = A*T^ Exp(-Ea/RT). Units are in moles, cubic centimeters, seconds, 


Kelvins, and calories/mole. 
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Table L Reaction Mechanism for Mg-0 2 


Reaction 


A 

P 

E a 

Ref. 

Mg + 02 = MgO + 0 

(Rl) 

4.44E+12 

0.5 

30500.0 

9 

Mg + O + M = MgO + M 

(R2) 

1.90E+14 

0.5 

0.0 

9 

MgO = MgO(l) 

(R3) 

6.11E+1 1 

-2.0 

17268.6 

10 

MgO(l) = MgO(s) 

(R4) 

l.OOE+15 

0.0 

0.0 

” 

20 + M = 0 2 + M 

(R5) 

1.89E+I3 

0.0 

-1788.0 

12 

The rate coefficients are in the form Kj ■ 
Kelvins, and calories/mole. 

Table 2. 

= A* 7^ Exp(-Ea/RT). Units are in moles, cubic centimeters. 
Reaction Mechanism for Mg-Air 

seconds, 

Reaction 


A 

P 

e a 

Ref. 

Mg + 0 2 = MgO + 0 

(Rl) 

4.44E+12 

0.5 

30500.0 

9 

Mg + O + M = MgO + M 

(R2) 

1.90E+14 

0.5 

0.0 

9 

MgO = MgO(l) 

(R3) 

6.11E+11 

-2.0 

17268.6 

10 

MgO(l) = MgO(s) 

(R4) 

l.OOE+15 

0.0 

0.0 

— 

20 + M = 0 2 + M 

(R5) 

1.89E+13 

0.0 

-1788.0 

12 

N + NO = N 2 + 0 

(R6) 

3.50E+13 

0.0 

330.0 

11 

N + 0 2 = NO + 0 

(R7) 

2.65E+12 

0.0 

6400.0 

11 

N 2 0 + 0 = n 2 + 0 2 

(R8) 

1.40E+12 

0.0 

10810.0 

11 

N 2 0 + O = 2NO 

(R9) 

2.90E+13 

0.0 

23150.0 

11 

N0 2 + 0 = NO + 0 2 

(RIO) 

3.90E+12 

0.0 

-240.0 

11 


The rate coefficients are in the form = A *7^ Exp(-Ea/RT). Units are in moles, cubic centimeters, seconds. 


Kelvins, and calories/mole. 






Table 3. Reaction Mechanism for Mg-C0 2 


Reaction 


A 

P 

e a 

Ref. 

Mg + C0 2 = MgO + CO 

(Rll) 

2.00E+14 

0.5 

34847.0 

1 

Mg + CO = MgO + C 

(R12) 

2.00E+1 1 

0.5 

30000.0 

1 

Mg + 0 2 = MgO + 0 

(Rl) 

4.44E+12 

0.0 

30500.0 

9 

Mg + O + M = MgO(s) + M 

(R2) 

1.90E+14 

0.0 

0.0 

9 

MgO = MgO(l) 

(R3) 

6.11E+11 

-2.0 

17268.6 

10 

MgO(l) = MgO(s) 

(R4) 

l.OOE+15 

0.0 

0.0 

— 

20 + M = 0 2 + M 

(R5) 

1.89E+13 

0.0 

-1788.0 

12 

co + o + m = co 2 +m 

(R13) 

6.02E+14 

0.0 

3000.0 

11 

CO + O 2 = CO 2 + 0 

(R14) 

2.50E+12 

0.0 

47800.0 

11 

c + o 2 = CO + 0 

(R15) 

5.80E+13 

0.0 

576.0 

11 


The rate coefficients are in the form K s = A*r Expi-E^RT). Units are in moles, cubic centimeters, seconds, 


Kelvins, and calories/mole. 



